Reactive Hall response 
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common expression used to describe the Hall constant 
in quasi-one dimensional conductors described by a band 
picture , (iv) good accord with Bh for ladder systems 
calculated using the numerical method proposed in 

The Hamiltonian In the following we will consider a 
generic Hamiltonian for fermions on a lattice, where for 
simplicity we describe the kinetic energy term by a one 
band tight binding model; it is straightforward to extend 
this formulation to a many-band or continuum system. 
The sites are labeled l{m) along the a;(?/)-direction with 
periodic boundary conditions in both directions: 
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The zero temperature Hall constant Rh, described by reactive (nondissipative) conductivities, is 
analyzed within linear response theory. It is found that in a certain limit, Rh is directly related 
to the density dependence of the Drude weight implying a simple picture for the change of sign of 
charge carriers in the vicinity of a Mott-Hubbard transition. This novel formulation is applied to 
the calculation of Rh in quasi-one dimensional and ladder prototype interacting electron systems. 

PACS numbers: 71.27.4-a, 71.10.Fd, 72.15.Gd 

It is now well known that in strongly correlated sys- 
tems the, zero temperature (T=0), reactive part of the 
conductivity can be used as a criterion of a metallic or 
insulating ground state Q. In particular, following the 
work of Kohn, the imaginary part of the conductivity, 
cr"{uj — > 0) = 2D/ijj, characterized by D (now called the 
"Drude weight" or charge stifness), can be related to the 
ground state energy density dependence on an applied 
fictitious flux as D = {l/2)d'^e^ /d(l)%^Q. 

A similar question is posed by the doping of an in- 
sulating state, where it would be interesting to have a 
simple description of the charge carriers sign as probed 
in a Hall experiment. For instance, we would like to de- 
scribe the doping of a Mott-Hubbard insulator; within a 
semiclassical approach it is expected that the Hall con- 
stant Rh — +l/e(5, hole-like (positive) near half-filling 
{5 = 1 — n, n=density), changing to Rh — — 1/en, 
electron-like at low densities, the turning point depend- 
ing on the interaction. 

Over the recent years, ingenious ways have been 
proposed |^,^ for characterizing this sign change and 
strongly correlated electron systems, as the t — J model 
have been studied. In particular, following the suggestion 
to focus at the T=0 Hall constant within linear response 
theory |^ , the Rh of a hole in the t — J model was ana- 
lyzed and a numerical method was proposed for calculat- 
ing the Hall response in ladder systems Q . This activity 
is partly motivated by the physics of high temperature 
superconductors viewed as doped Mott-Hubbard insula- 
tors and related Hall measurements showing a change of 
the sign of carriers with doping . 

In this work, we will show that within a certain 
frequency-w, wavevector-g limiting procedure, the T=0, 
w — *■ 0, thus "reactive" Hall constant, is simply related to 
the density dependence of the Drude weight. Following 
this point of view, we recover in a straightforward way: 
(i) the semiclassical expressions for Rh at low density and 
near an insulating state, (ii) a physical picture of the sign 
change of carriers in the vicinity of a Mott-Hubbard tran- 
sition and its dependence on interaction strength, (iii) a 
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Q,Tn(c|„) is an annihilation (creation) operator at site 
(l,m) and the spin is neglected as it enters in a trivial 
way in the formulation. The U term can represent a 
many-particle interaction or a one particle potential. We 
take the lattice constant so as to consider a unit volume, 
electric charge e = 1 and h = 1. We add a magnetic 
field along the z-direction, modulated by a one compo- 
nent wavevector-g along the y-direction, generated by the 
vector potential Am ; this allows to take the zero magnetic 
field limit smoothly B: 
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(for convenience, we will present the long wavelength 
limit, substituting 2sin(g/2) q). Electric fields along 
the x,y directions are generated by time dependent vec- 
tor potentials: 
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£;^(t) = i;^e-*^*, Ey{t)=iEye-'^^- z = Lo + ir]. (3) 

Currents are defined through derivatives of the Hamilto- 
nian expanded to second order in 0^'^: 

with the paramagnetic parts: 

l,m 

jy ^t'J2 e^^("+i/^) (*ct ™+iQ,™ + h.c). (5) 

The reactive Hall response From standard hnear re- 
sponse theory we obtain: 

{J-)^a,.,.E%t)+a^.^yEyit) 
{jy)^a,y,.Em)+a^y^.Ey{t) . (6) 

(...) are ground state expectation values in the presence 
of the magnetic field, with the conductivities 
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Now, in contrast to the usual derivation of the Hall con- 
stant expression, we will keep the dependence explicit 
by converting the current-current to current-density cor- 
relations using the continuity equation: 

{JD - --Xn,,^E^{t) + C-)\ '-Ey{t) , (8) 
H q q ^ ^ z 

with Uq = Ei,m(-*e''^'")4m^'.'"- 

At T=0, the response is non-dissipative so we will 
study the reactive (out-of phase) induced currents. Fur- 
thermore, at this point we will consider the "screening" 
(or slow) response in the y— direction, by taking the {q, lo) 
limits in the order u —>■ first and q — > last; in the usual 
"transport" (or fast) response the limits are in the oppo- 
site order ||^. As we will discuss below, this approach 
leads to a simple physical picture for the Hall constant 
and it might be argued that at least for certain cases, 
for example for a system of finite size in the direction, 
it is indeed the right one. The expressions (^ for the 
currents become: 

(J-)o = <.,.(^-0))(*i?-(i)) 

+ -x;-„,(^-o)i?^(t) 

q ■< 1 

+ {-r-Xn,^^u;^0){^Ey{t)) , (9) 
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where the subscript zero denotes the leading order in uj 
response. 
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and \n){En) are eigenstates (eigenvalues) of the Hamil- 
tonian in the presence of the magnetic field. 

Now, following Kohn's observation |^ , we can identify 
the different terms as derivatives of the ground state en- 
ergy density e° of a fictitious Hamiltonian depending on 
static (p^ , fjiq fields: 

H = i-t) + h.c.) 

l,m 

+ i-t') + h.c.) + ^iqnq + U. (11) 
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For iJ(A,/i), using the following identity, 
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we can rewrite the currents as 
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{Jlo = ^-^{^E4t)) + iJ-)e%Ey{t) 



iJ> = -q^U^E4t)~-el^S'Ey{t)). (13) 

Finally, setting (J|)o = we determine the "reactive" 
Hall constant: 
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Neglecting the cross-terms e° ^xC^^;^ and Taylor ex- 
panding the numerator in _B, we can rewrite Rh as: 
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Using (|T2|) we find the final expression: 



Rh — — 



Dk„ 



(15) 



(16) 



where, 
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= i-t) T.(-'^''"')(^ll,r.^^'^^"^ + h.c). (17) 
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the Drude weight, is identical to Dq by 
the replacement of (T^) by f (T^). = ej^_^^^ = 
duq/dfiq is the compressibility corresponding to the den- 
sity modulation Uq. Notice that the spatial dependence 
of jq and Uq is the same as that of Am ■ 

Taking the q ^ limit, we obtain a particularly simple 
expression for Rh- 
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A handwaving argument leading to expression ( p^ for 
i' ^ is as follows: Am corresponds to a twist of bound- 
ary conditions on chain— to, inducing an extra current on 
each chain proportional to D (besides the uniform one in- 
duced by the flux (j)^); minimization of the energy at fixed 
X— current gives rise to an to— dependent charge density. 
This induced charge density can then be canceled by the 
"Hall potential" fiq |^ . Note that a similar idea, analyz- 
ing the Hall constant in terms of independent channels 
(edge states) , exists in the literature of the Quantum Hall 
effect 0. 

This expression is appealing as it gives a direct, intu- 
itive understanding for the change of sign of charge car- 
riers in the vicinity of a metal- insulator transition. First, 
at low densities, D n giving Rh — —\/n] close to a 
Mott insulator D oa 5 — 1 — implying Rh — +1/(5. 
Furthermore, we obtain a change of sign in the vicin- 
ity of a Mott transition at a density which depends on 
the interaction strength and is given by the position of 
the maximum of D. Second, for independent electrons, 
where D is proportional to the kinetic energy, by taking 
the limit i' — » and calculating I? as a sum of I?'s for 
individual cc— chains, we obtain from (fi^): 
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an expression used for the Hall constant of quasi one- 
dimensional compounds . Considering that the t' — > 
limit might by subtle, it is of particular theoretical and 
experimental interest whether the Hall constant of quasi- 



one dimensional correlated systems |11 is indeed given 
by the expression and thus related to the Drude weight of 
the individual chains. The same applies for the transverse 
Hall effect of weakly coupled planes. 

Examples In this section we present a generic pic- 
ture for the behavior of the Hall constant for models 
of strongly correlated fermions showing a Mott-Hubbard 
metal- insulator transition. This picture emerges, on the 
one hand, by an exact calculation of Rh for ladder sys- 
tems using the numerical method of ref. |^ and on the 
other hand, from the expression (^|) assuming nearly 
decoupled chains [t' 0) and calculating D{n) for each 
chain analytically using the Bethe ansatz method ||l2|,|l3| . 
It is clear that this analytical approach refers to either 
ladder (with t' 0) or quasi one- dimensional models. 



Three prototype models will be discussed: the Hub- 
bard model, as the most experimentally relevant, the 
spinless fermions model ("t-V") showing both a metal- 
lic and an insulating phase depending on interaction 
strength and the supersymmetric t — J model. 

(i) The Hubbard model is given by the Hamiltonian: 

H= i-t)J2(^Ul,rn,aCl,m,a+h.C.) 

l,m 

(20) 
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ci,m,cr{cl m ij) annihilation (creation) operator at site 
(I, to) of a fermion with spin a =]" , |. Rh extracted from 
a Bethe ansatz calculation of D{n) for the one dimen- 
sional Hubbard model is shown in Fig. 1. 




FIG. 1. Rh for the Hubbard model from expression 
for t' 0. 

This behavior is characteristic of correlated systems 
undergoing a metal-insulator transition at half-filling: 
at low densities Rh — ^^/n, while near half- filling 
Rh — +1/^, the position of change of sign of the car- 
riers depending on the details of the interaction. 

(ii) The t- V model on a ladder is given by: 



l.m 

(21) 

I l,m 

Here and in the following / = 1, Lx,m — 1,2. For a 
single chain, this model describes a metallic phase at all 
densities for V < 2t, while for V > 2t it is an insulator at 
half-filling. In Fig. 2 we show Rh calculated numerically 
on finite systems for two values of t' and analytically from 
( p8| ) in the t' — > limit. The numerical evaluation being 
especially sensitive to finite size effects for t' — > 0, we 
study relatively large values of t'. 

Results for Rh clearly show the difference between the 
metallic regime V = t, where at half-filling {n — 0.5) we 
get Rh — 0, while in the insulating regime V — At, we 
are dealing with Rnin — > 0.5) oo. 
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FIG. 2. i?H for the t — V ladder from expression (|18|) for 
— > (continuous line) and from a numerical evaluation 
(symbols). V — t{4:t), metallic (insulating) phase ai n — 0.5. 

(iii) The t-J model on a ladder is given by the Hamil- 
tonian: 



H = {-t)J2^'^Ul,7n.aCl,^^,<r + h.C.) 
+ {-t')Y.^4x-CU2,a + h.C.) 



l.m 
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Sim is the spin operator at site (/, m) and the double 
occupancy on a site is forbidden. 




FIG. 3. Rh for the t — J ladder from expression (g^ for 
t' ^ Q (continuous line) and from a numerical evaluation 
(symbols) . 



the "supersymmetric" model, J = 2t, and by numerical 
evaluation for t' = 0.5t and different size systems. 

The above three examples show a remarkable agree- 
ment between the numerical evaluation of on finite 
size systems using the method of ref. |^ (at finite t') and 
the analytical calculation using (|l|) for i' 0, indicat- 
ing a relative insensitivity on the transverse coupling t' 
for ladders. These results confirm the intuitive picture 
for the behavior of the Hall constant in the vicinity of a 
metal-insulator transition and present an intriguing link 
between the Hall constant and the Drude weight. It is 
possible that Rh is dominated at low temperatures by 
correlations and not the relaxation mechanism so this 
formulation could have more general validity. 

In conclusion, the emerging simple physical picture 
raises the question of the relation of this novel formula- 
tion to the traditional semiclassical approach to the Hall 
constant, its range of validity, the role of relaxation in 
the description of the Hall effect and of the perspectives 
for an extension at finite temperatures. 
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In Fig. 3 we show again Rh calculated analytically for 
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